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Topological phases of matters are of fundamental interest and have promising applications. Fasci¬ 
nating topological properties of light have been unveiled in classical optical materials. However, the 
manifestation of topological physics in quantum optics has not been discovered. Here we study the 
topological phases in a two-dimensional momentum-space superradiance lattice composed of timed 
Dicke states (TDS) in electromagneticahy induced transparency (FIT). By periodically modulating 
the three EIT coupling helds, we can create a Haldane model with in-situ tunable topological prop¬ 
erties, which manifest themselves in the contrast between diffraction signals emitted by superradiant 
TDS. The topological superradiance lattices provide a controllable platform for simulating exotic 
phenomena in condensed matter physics and offer a basis of topological quantum optics and novel 
photonic devices. 


The quantum Hall effect (QHE) [T] reveals a new class of 
matter phases, topological insulators (TIs), which have 
been extensively studied in solid-state materials Hd 
and recently in photonic structures IHHI2, time-periodic 
systems HMS] and optical lattices of cold atoms HZ]. 
The first TI is the Haldane model [2] proposed in 1988, 
which shows that the QHE is an intrinsic topological 
property of the energy bands. It inspired the discov¬ 
ery of the quantum spin Hall effect mg and topo¬ 
logical superconductors [7]. The Haldane model con¬ 
sists of a honeycomb tight-binding lattice with complex 
next-nearest-neighbour (NNN) hopping, which breaks 
the time-reversal symmetry and induces an energy gap 
between two bands that have opposite Chern numbers. 
Inversion symmetry breaking by an on-site potential can 
also open band gaps. The interplay of these two types of 
symmetry breaking leads to transitions between phases 
with Chern numbers 0 and ±1. 

Notwithstanding its foundational role in topological 
condensed matter physics, the Haldane model has never 
been realized in solid-state systems. Floquet modulation 
of a circular polarised light in graphene was proposed to 
induce the complex NNN tunnelling HSHig. However, 
the required light is soft X-ray which would directly ex¬ 
cite electrons and be absorbed. The Haldane model of 
cold atoms in optical lattices m were recently realised 
in experiments. 

Our recent study shows that Scully’s timed Dicke 
states (TDS) [18] can form a superradiance lattice (SL) 
in momentum space m- Here we propose a quantum 
optics realisation of an in-situ tunable Haldane model us¬ 
ing two-dimensional SLs in a simple electromagneticahy 
induced transparency (EIT) configuration [20|. Since 
Dicke’s seminal paper in 1954 [21], superradiance has 




Figure 1 : Realisation of the Haldane model in a superradi¬ 
ance lattice of timed Dicke States in electromagneticahy in¬ 
duced transparency, (a) Schematic configuration of the probe 
field kp, EIT coupling fields along ki/2/3 and diffraction fields 
along k+ = kp + ki — ka and k_ = kp + ki — k2. (b) The en¬ 
ergy level diagram of the EIT coupling, probe and scattering 
fields, (c) Honeycomb structure of the SL. |ek) and |mk) cor¬ 
respond to the two sublattices. The NNN hopping D31 (red 
solid arrow), D23 (blue dot arrow) and D12 (green dash dot 
arrow) enclose a nonzero effective magnetic hux in momen¬ 
tum space. The arrow on the circle denotes the direction of 
excitation current when the Chern number C = 1 . 


been an important topic in quantum optics. A single 
photon with wave vector kp can excite a spatially ex¬ 
tended A/'-atom ensemble from the ground state \G) = 
151,52, • • • ,5iv) to the TDS [E] 


N 

ekp) = (1) 
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Figure 2 : Topological phases of the superradiance lattice, 
(a) The value of a (defined in Eq. if) as a function of the 
modulation phases of the coupling fields 02 and 03. 0 i = 0 . 
The signs “±” are for the Chern numbers C = ± 1 . The 
dashed line separates the two topological phases, (b) The two 
Chern numbers corresponding to two topologically distinctive 
configurations of the three phase angles 0 i, 02 and 03 on a 
unit circle. 


where and gi are the excited and ground states of the 
ith atom at position r^, respectively. The TDS stores 
a light momentum Jikp via phase correlations between 
atoms excited at different positions. This momentum 
can be transferred back to a single photon via directional 
emission [18]. By coupling |e) to another state \m) with 
three coherent plane wave fields, we construct a honey¬ 
comb SL of TDS in momentum space m, as shown in 
FigO The SL Hamiltonian with rotating-wave approxi¬ 
mation is (see Supplementary Information) 

H =—^ - |mk-ki)(mk-ki|) 

k 

3 (^) 

-EE fifl/[|ek)(mk-kj + h.c.], 

k 1=1 

where k = kp + r(k 2 — ki) + ^(ks — k 2 ) with integers 
r and s, fti is the Rabi frequency of the coupling field 
along k^, |mk) are defined the same as in Eq. 0 but for 
\m) states, and Ac = ujem — is the detuning of the 
transition frequency cjem between |e) and \m) from the 
angular frequency Uc of the EIT coupling fields. The 
states |ek) and |mk-ki) form two sublattices of a honey¬ 
comb structure in momentum space [see Fig. Be)]. 
causes the nearest-neighbour hopping. As we will show 
below, periodically modulating Qi can introduce NNN 
hopping with controllable phases. On the other hand, 
the detuning Ac causes an energy offset between the two 
sublattices, breaking the inversion symmetry. Thus three 
Haldane phases with Chern numbers ±1 and 0 can be in- 
situ realised. 

We first show how to induce complex NNN hopping. 


We consider periodically modulated Rabi frequencies Qi 
(/ = 1,2,3) of coupling fields, 

fli =^s + 2f2d cos{vdt -^4>i), (3) 


where and Q.d are the static and dynamic compo¬ 
nents of the Rabi frequency, Vd is the modulation fre¬ 
quency and 0/ is the modulation phase. Here we have 
assumed the size of the atomic ensemble be much smaller 
than cjvd (where c is the speed of light) and hence ne¬ 
glected the position dependence of the modulation. On 
the other hand, the ensemble size is much larger than 
cjvc and the number of atoms A ^ 1 such that the 
TDS in the SL are approximately orthogonal to each 
other, i.e., (ek'|ek) ~ ^kk'- The periodic modulation in¬ 
duces Floquet sidebands [2] with energy separation 
While Os is the intra-sideband nearest-neighbour hop¬ 
ping, fid induces the inter-sideband hopping. We choose 
Ud ^ so that the Floquet sidebands are well sep¬ 

arated in energy. Thus by second-order perturbation, the 
effective Hamiltonian of the NNN intra-sideband transi¬ 
tion mediated by Cld is (see Supplementary Information) 
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H = ^0/j(|ek+ki-kj)(6k| 

k i^j=i w 

T Inik—ki ) (nik—ki+k; —kj I): 

where the NNN hopping coefficient Q^ij = 2iO' sin(0/—0^) 
with O' = The fact that the NNN hopping co¬ 

efficient is purely imaginary is crucial to the topological 
phases. The loop transitions via NNN hopping accumu¬ 
late nonzero phases, as shown in Fig{^ (c). H' opens a 
band gap e^r = 4\/3fi|(a|0' where a is a dimensionless 
quantity determined by the summation of 

a =-Esin(</'i+i-</>;) =4nsinA2_^. (5) 

1 = 1 1 = 1 


The Chern numbers of the upper and lower bands C and 
C' are C = —C' = sign((a) (see Supplementary Infor¬ 
mation). In Fig. |S.2| (a), we plot a with 0i = 0 and 
0 < 02,3 < 27r. The topological property of this SL Hal¬ 
dane model can be represented by the distribution of 0/ 
on a unit circle. There are two distinct topological con¬ 
figurations, counter-clockwise 0i, 02 and 03 for (7 = —I, 
and clockwise 0i, 02 and 03 for C = I, as shown in 
Fig |S.2| (b). The time reversal t ^ —t in Eq.(S.4) is 
equivalent to 0/ ^ —0/, which leads to C —C. 

Unlike TIs [5j |6], topological superradiance lattices 
(TSLs) have no outer edges in the semiclassical limit 
of the coupling fields (see Supplementary Information). 
Neither do TSLs have Fermi surfaces. Nonetheless, the 
TSL has its unique topological properties that are observ¬ 
able. The TDS have directional superradiance emission. 
Of all the TDS in the SL, only those |ek) with c|k| — ujeg 
within the energy bands {uJeg is the transition frequency 
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Figure 3: Topological phase transitions signatured by super¬ 
radiance contrast, (a) The contrast 77 (defined in equation (6)) 
as a function of the modulation phases 02 and 03 with 0i = 0. 
The red dot lines are zero points. Ap = Ac = 0. Q' = 0.01. 
(b) The topological phase transition between C = 1 and 0. 
Ap = Ac/2. Q' = 0.01 (red), 0.1 (green), 0.5 (blue) and 1.0 
(black), (c) The SL band structures in the topological phase 
transition. = 0.1. K and K' are the two valley points. For 
(b) and (c), 0i = 0, 02 = 47r/3 and 03 = 27r/3. For all three 
figures, = 3, 7e = 1 and 7^ = 0.1. 


between |e) and \g)) can satisfy both energy and mo¬ 
mentum conservation, and have directional emission in k 
m- We call these states superradiant TDS and the other 
ones subradiant TDS. We can regard these superradiant 
TDS as an inner edge of a honeycomb lattice of subra¬ 
diant TDS [23]. The topological orders lead to different 
light emissions from different superradiant TDS. Alterna¬ 
tively, we can also tune the probe field frequency to test 
the topological band properties at certain energy, which 




Figure 4: Real-space propagation of a light pulse through 
an atom ensemble driven into a trivial or topological super¬ 
radiance lattices. Diffraction of a weak Gaussian probe pulse 
in x-y plane for (a) a trivial SL with = 0 and (b) a 
TSL with C = 1. 01 = 0, 02 = 47r/3, 03 = 27r/3 and 

Q' — 0.01. The derived linear and nonlinear susceptibili¬ 
ties are — ^0.1410, x+ = X- = —z0.0522 in (a) and 
X^^^ = z0.1446, x+ = —z0.1057, x- = —i0.0021 in (b). The 
parameters are such that Qg = 3, Ap = Ac = 0, 7e = 1 
and 7m = 0.1. The square denotes the area occupied by 
atoms. The three figures in each group are for scaled time 
t = 1,80,100 successively. We assume the group velocity be 
uniform everywhere. 


is analogue to tuning the Fermi surface in a fermionic 
system. 

For the sake of simplicity, we set the wavevectors of 
the three EIT coupling fields to be ki = —kcX^ k 2 = 
kc{x — V3y)l2 and ks = kc{x \/3^)/2, and the probe 
field wavevector k^ = —ki. In this case we have only 
three superradiant TDS with wavevectors kp, k+ = — ks 
and k_ = — k 2 , and the diffraction fields are along k±, 
as shown in Fig0 (a). We set all fields on resonance, i.e., 
Ac = 0 and Ap = ujeg — i^p = 0. The excitation |ekp) 
flows to |ek±) and emits photons along k±. We denote 
the steady state probability amplitudes of states |ek±) as 
Ck± and define the superradiance contrast 


CkJ^-|Ck_| 

Ck+P + |ck_ p 


For C = 1, the excitation current flows along jek^) ^ 
|ek+) ^ Pk_) [23], as shown in Figj^ (c). Since each 
TDS pk) or |mk) has decoherence rate 7e or 7 ^, respec¬ 
tively, the excitation decays while flowing and it is more 
probable in state pk^) than in state pk_)- We therefore 
have 77 > 0. Similarly, for C = — 1, 77 < 0. Thus the sign 
change of the superradiance contrast signatures the topo¬ 
logical phase transition, as seen in Fig|3|(a). The super¬ 


(b) except for the two diagonal corners where the topo¬ 
logical currents are weak and the local effect inside a unit 
cell dominates (see Supplementary Information). 

The topological phase with C = 0 for non-zero can 


radiance contrast in FigJ^ (a) is consistent with Fig S.2 
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be reached by breaking the inversion symmetry. Sub¬ 
stantially easier than in graphene electrons, a sublattice 
offset in SLs can be introduced by choosing Ac 7 ^ 0. We 
set (pi = 47 r(/ — l)/3 and thus C = 1 when Ac = 0. The 
energy gap in absence of sublattice offset is e^r = IShQ'. 
Topological phase transitions occur at h\Ac\ = etr [ 2 ] (see 
Supplementary Information). In FigJ^(b), we plot 77 as 
a function of hAc/etr- For = 0.01, the phase transi¬ 
tion is smeared out by the relatively large 7 e. The phase 
transition gets more apparent for larger Q'. Depending 
on 77 has peaks or kinks near the phase transition 
point because the two bands touch at the middle of the 
band gap, where the probe field probes, as shown in Fig|^ 
(c). 

The TSLs have unique features in transient light prop¬ 
agation under pulse probe. In Figj^ we compare the 
pulse propagation in a trivial SL with zero ft' and in a 
TSL. For a weak probe pulse, the linear susceptibility is 
oc Ck and the linear absorption is Imy^^^ The two 
multi-wave-mixing signals along k± correspond to the 
nonlinear susceptibilities oc Ck± and can be under¬ 
stood as a result of optical grating m- We simulate the 
pulse propagation for the three modes along using 
coupled wave equations [24] (see Supplementary Infor¬ 
mation). For a trivial SL without modulation, the light 
propagating along k± is symmetric, while for a TSL with 
C = 1 , the topological currents drive the probe pulse to 
k+, even if the NNN hopping is two orders of magnitude 
smaller than the nearest-neighbour hopping. 

TSLs can be readily realised in experiments for cold 
alkali atoms. Taking ^^Rb D1 line for example, we can 
have l^f) = |5^5 'i/2,T = 2 ), |e) = |5^Pi/2,T = 2 ) and 
\m) = | 5 ^ 5 'i/ 2 ,P = 3). 7 e = 2.9MHz anci 7 ^ is control¬ 
lable via inhomogeneous magnetic field [25|. The Rabi 
frequency = 37e = 8 . 6 MHz (intensity 25mW/cm^). 
The modulation frequency can be Ud = lO^e = 28.8MHz 
which is large enough to separate the Floquet bands. One 
can trap 10 ^ atoms in 1 mm^ such that A ^ 1 and 
the size L of the ensemble cjv^ <C L cjvd- fhe 
jiK regime, the thermal random motions have negligible 
Doppler shifts (^kHz). Another possible type of phys¬ 
ical systems are rare earth atoms doped in solids [26] . 
One should first optically pump nearly all population to 
1^), then turn on three optical fields coupling |e) to |77i), 
send in a weak field probing the \g) to |e) transition, and 
detect the diffraction signals. 

Similar to the electronic TIs [musi and the optical 
lattice simulations of TIs the topological properties 
of TSLs are determined by Schrodinger equations, dif¬ 
ferent from the photonic TIs isHm, which are governed 
by Maxwell equations. The generation and detection of 
the topological properties of TSLs, however, can be eas¬ 
ily controlled by light. The unique feature of the TSL is 
that its lattice sites have discrete momenta rather than 
positions. It has the advantage to be extended to dimen¬ 
sions higher than three where no real space lattices exist 


m and offers a platform for high-dimensional topologi¬ 
cal physics m- 


Funding Information 

Science Foundation Grants No. PHY-1241032 (IN¬ 
SPIRE CREATIV) and PHY-1068554; Robert A. Welch 
Eoundation (Grant No. A-1261); Herman E. Heep and 
Minnie Belle Heep Texas A&M University Endowed 
Eund; Hong Kong RGG/GRE; GUHK VG’s One-Off 
Discretionary Eund; Natural Science Eoundation of 
Ghina U1330203. 

The authors thank Marian O. Scully for helpful dis¬ 
cussions. 


Electronic address: cuhkwdw@gmail.com 
^ Electronic address: rbliu@phy.cuhk.edu.hk 

[1] K. von Klitzing, G. Dorda and M. Pepper, Physical Review 
Letters 45, 494-497 (1980). 

[2] F. D. M. Haldane, Physical Review Letters 61, 2015-2018 
(1988). 

[3] C. L. Kane and E. J. Mele, Physical Review Letters 95, 
226801 (2005). 

[4] B. A. Bernevig and S. C. Zhang, Physical Review Letters 
96, 106802 (2006). 

[5] M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buhmann, 
L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 
318, 766-770 (2007). 

[6] C.-Z. Chang, J. Zhang, X. Feng, J. Shen, Z. Zhang, M. 
Guo, K. Li, Y. On, P. Wei, L.-L. Wang, J.-Q. Ji, Y. Feng, 
S. Ji, X. Chen, J. Jia, X. Dai, Z. Fang, S.-C. Zhang, K. 
He, Y. Wang, L. Lu, X.-C. Ma, Q.-K. Xue, Science 340, 
167-170 (2013). 

[7] X.-L. Qi, T. L. Hughes, S. Raghu, S.-C. Zhang, Physical 
Review Letters 102, 187001 (2009). 

[8] F. D. M. Haldane and S. Raghu, Physical Review Letters 
100, 013904 (2008). 

[9] M. Hafezi, S. Mittal, J. Fan, A. Migdall and J. M. Taylor, 
Nature Photonics 7, 1001-1005 (2013). 

[10] M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, D. 
Podolsky, F. Dreisow, S. Nolle, M. Segev and A. Szameit, 
Nature 496, 196-200 (2013). 

[11] A. B. Khanikaev, M. S. Hossein, W.-K. Tse, M. Kargar- 
ian, A. H. MacDonald and G. Shvets, Nature Materials 
12, 233-239 (2013). 

[12] A. V. Nalitov, G. Malpuech, H. Tergas and D. D. Sol- 
nyshkov. Physical Review Letters 114, 026803 (2015). 

[13] T. Oka and H. Aoki, Physical Review B 79, 081406 
(2009). 

[14] J.-I. Inoue and A. Tanaka, Physical Review Letters 105, 
017401 (2010). 

[15] T. Kitagawa, E. Berg, M. Rudner and E. Dernier, Phys¬ 
ical Review B 82 , 235114 (2010). 

[16] N. H. Lindner, G. Refael, and V. Galitski, Nature Physics 
7, 490-495 (2011). 



5 


[17] G. Jotzu, M. Messer, R. Desbuquois, M. Lebrat, T. 
Uehlinger, D. Greif and T. Esslinger, Nature 515, 237-240 
(2014). 

[18] M. O. Scully, E. S. Ery, G. H. R. Ooi and K. Wodkiewicz, 
Physical Review Letters 96, 010501 (2006). 

[19] D.-W. Wang, R.-B. Liu, S.-Y. Zhu and M. O. Scully, 
Physical Review Letters 114, 043602 (2015). 

[20] K.-J. Boiler, A. Imamoglu and S. E. Harris, Physical 
Review Letters 66, 2593-2596 (1991). 

[21] R. H. Dicke, Physical Review 93, 99-110 (1954). 

[22] J. H. Shirley, Physical Review 138, B979-B987 (1965). 

[23] Y. Lumer, Y. Plotnik, M. G. Rechtsman, and M. Segev, 
Physical Review Letters 111, 243905 (2013). 

[24] D.-W. Wang, H.-T. Zhou, M.-J. Guo, J.-X. Zhang, J. 
Evers and S.-Y. Zhu, Physical Review Letters 110, 093901 
(2013). 

[25] V. B. Tiwari, S. Singh, H. S. Rawat, M. P. Singh and S. 
G. Mehendale, Journal of Physics B 43, 095503 (2010). 

[26] G. W. Thiel, T. Bottger, and R. L. Gone, Journal of 
Luminescence 131, 353-361 (2011). 

[27] S.-G. Zhang and J. Hu, Science 294, 823-828 (2001). 



Figure S.l: Timed Dicke states and the associated photon 
numbers in an SL. The numbers associated with each TDS 
are the photon number difference from those of state |ekp). 


Supplementary I 


EFFECTIVE HAMILTONIAN 


The three-level atoms we use to construct the two- 
dimensional (2D) superradiance lattice (SL) have a 
ground state \g)^ an excited state |e) and a third state 
\m). The optical fields that couple |e) and \m) have three 
modes with wave vectors ki = —kcX^ k 2 = kc{^x — ^y) 
and ks = kc{^x ^y). The interaction Hamiltonian 
with rotating-wave approximation is 

3 N 

He = + H.c., (S.l) 

1=1 i=i 

where ki and ai are the vacuum coupling strength and 
annihilation operator of the Ith mode, respectively. N 
is the number of atoms and vj is the position of the jth 
atom. 

Initially the atomic ensemble is in the ground state 
\G) = l^fi,^ 2 , •••,^at)- a weak probe field with wave vec¬ 
tor kp can prepare the atomic ensemble in the timed 
Dicke state (TDS) 




1 ^ 

) = 7f (S.2) 


i=i 


by the collective absorption of a single photon. The 
combinational quantum states of the atoms and three 
coupling modes can be written as |ekp,ni,n 2 ,ns), where 
ni (/ = 1,2,3) is the photon number in mode k/. This 
state is coupled to |mk -k 2 5 '^2 + 1, ns) with coupling 

strength —fi/^ 2\/^2 + 1- The state |mk -ks) is defined 
the same as in equation (S.2) with Cj replaced by ruj 


and kp replaced by k^ — k 2 . Similar coupling also ex¬ 
ists for the other two modes. Inikp-ks ?'^2 + l,ns) 
is in turn coupled to jck^+ki-ks?~ 1,^2 + l,ns) or 
|ekp+k 3 -k 2 5 ^2 + l,ns — 1) via excitation by modes 

ki or ks- These states form a honeycomb lattice with 


discrete momentum coordinates, as shown in Fig. S.l 


called the superradiance lattice [T]. The two sublattices 
of the SL correspond to TDS for |e) and |m). We de¬ 
note the coupling field frequency as Uc and the transition 
frequency between |e) and |m) as co’em- The energy dif¬ 
ference between the two sublattice is Ac = oJem — We 
set the zero energy at the middle of the energies of the 
two sublattices. Then the energies of the |e) and |m) 
sublattices are Ac/2 and — Ac/2, respectively. 


Although the coupling strengths in the SL are site- 
dependent, if we use coherent coupling fields with large 
average photon numbers {ni) ^ 1, the Rabi frequency of 
the Ith mode can be approximated as the classical Rabi 
frequency 0/ = ^z\/(^z)- The SL with quantized photon 
numbers has edges when either of the three coupling field 
photon numbers reduces to zero. The total lattice has a 
triangular boundary. However, since the TDS have finite 
life time, which we suppose to be r on average, once we 
create an excitation |ekp) in the SL, the average distance 
this excitation can travel is in the order of Qir. In this 
paper we have ftir {ni) and the edges can never be 
reached. Interesting edge effects might exist in the few- 
photon limit and it will be discussed elsewhere. Since 
the photon numbers are correlated to the momenta k’s 
in the TDS, we can simplify the notations by dropping 
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the photon numbers. The Hamiltonian of these TDS is 

^=^E(l^k)(ek| - |mk-k,)(mk-kj) 

3 (S.3) 

-EE fif^/[|ek)(mk-kzl + h.c], 

k 1=1 

where k = kp + r(k 2 — ki) + ^(ks — k 2 ) with integers r 
and s. 

The next-nearest-neighbour (NNN) hopping terms are 
induced by the periodic modulation of the three coupling 
field Rabi frequencies, 

9.i = Q.s + 29.d,cos(ydt-'k}^-v + 4>i), (S.4) 

where and Vtd are the static and dynamic components 
of the Rabi frequencies. Vd is the modulation frequency. 
(^i is the modulation phase of the Ith field, kj^ is the 
modulation wavevector. We choose the atomic ensemble 
much smaller than l/|k^| such that the position depen¬ 
dence of the phase —k[^ • r can be neglected. We set 
—k[^ • r = 0 for simplicity. On the other hand, the size of 
the atomic ensemble shall be much larger than Xjvc and 
the total number of atoms ^ 1 such that 

(ek'|ek) = ^Ee*(‘^-‘^')-^«4k'. (S.5) 

i=i 

We assume the atoms be randomly distributed and their 
number is large enough to cover all possible points in the 
real space Brillouin Zone. The SL can be regarded as 
infinite. 

We expand the Hamiltonian into static, positive- and 
negative-frequency components, 

H = + + (S.6) 

where 

Ho =^E(|ek)(ek| - |mk-k.)(mk-kj) 

3 (S-7) 

- Hfls EE[l®k)(”^k-k,| + /1-C.], 

k 1 = 1 

3 

= —hCid EE^*^' [|ek)(mk-k,|+/1-C.], (S.8) 

k 1=1 


The dynamics of the system is a Floquet problem EE]. 
Based on the Floquet theorem, the wave function can be 
written as 


1^-) = e--‘/'^|^(i)), 
where e is the quasi-eigenenergy and 


(S.IO) 


(S.ll) 


is a periodic wave function with period 27r/z/^. Substitut¬ 
ing equations (|S.6[)-(S.ll ) to the Schrodinger equation. 


ihd\'^)/dt = i7|T), we obtain 


n 

n 

(S.12) 

The terms with the same time evolution phase factors 
should be equal on both sides. We therefore have 


(e - hnVd)\i’n) = Ho\'tpn) + H+i\'tpn-l) + -H'-ilV’n+l)- 

(S.13) 

The quasi-eigenenergy can be obtained by diagonalizing 
the above Hamiltonian. 

For the sake of simplicity, we assume the separation 
between the Floquet sidebands is much larger than the 
bandwidth, Vd ^ where the perturbation the¬ 

ory can be applied ims]. When the probe field is near 
resonance, Ap = uj^g — Up fig where ooeg is the tran¬ 
sition frequency between |e) and \g) and Up is the probe 
field frequency, only the Floquet band with n = 0 in 


Eq.(S.13) is relevant. For states with eigenfrequencies 


near e = 0, the effective Hamiltonian can be obtained by 
standard second-order perturbation as 


elV’o) = Heslipo), 


(S.14) 


where 


Feff = Ho + H', 


(S.15) 


with NNN Hamiltonian 


(S.16) 


Explicitly, 


3 

= -hna E [|ek)(mk-k, I + h.c.]. (S.9) 

k 1=1 

Note that H±i are not Hermitian themselves, but H’+i 
is the Hermit ian conjugate of H_i. The phase factors 
H±i play the crucial role for the complex NNN 
hopping terms in the Haldane model. 


k lj^j=l 

(l^k+k;—kj)( 6 k| + I'^k—ki ) (^k—ki+ki—kj I) 5 

(S.17) 

where Vtij = = 2iO'sin(0/ — (j)j) 

with = O^^/ud- The crucial factor i = \/^ comes 
from the quantum interference between the two pathways 




7 
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(a) 


(b) 


Figure S.2: Complex next-nearest-neighbor hopping in¬ 

duced by Rabi frequency modulation, (a) The Floquet side¬ 
bands of the SL and the NNN transition Q 13 (dashed arrow) 
induced by the interference between two second-order inter¬ 
sideband transitions (solid arrow). Phases are labelled for 
each step, (b) The effective NNN transitions in a unit cell. 
The NNN hopping Q 31 (red arrow), Q 23 (blue dot arrow) and 
Qi 2 (green dash dot arrow) enclose a nonzero effective mag¬ 
netic flux. 


shown in Fig |S.2| (a). The loop transitions via NNN hop¬ 
ping accumulate nonzero phases due to these complex 
factors, as shown in Fig ]S.2| (b). 

The effective Hamiltonian is greatly simplified in the 
real-space representation. We denote the real-space basis 
states as 


\^rj) 1 ^ 1 : • • • 5 • • • 5 S'At), (S 18 ) 

l^r,-) = 1^1,^2, . .. . . . ,^Ar), 

The effective Hamiltonian can be written as 

= (S.19) 

3 

where the effective magnetic field h(rj) = {hx^hy^hz) 
with 



Figure S.3: Effective magnetic fields in a Brillouin zone 
for topological and trivial SL’s. (a) A topological SL with 
Ac = 0. (b) A trivial SL with Ac = 36QL The arrows show 
the vector hxX + hyy and the colours show hz- The x and y 
axis are in unit of k~^. The phase of the coupling field is 
(j)i = {I — l) 47 r /3 and a = 3\/3/2. Note that the KjK' points 
at the boundaries of the first Brillouin zone are equivalent. 


BERRY CONNECTION AND BERRY 
CURVATURE 


3 

^ cos(rj ■ ki), 
1=1 


(S.20) 


The topological properties of the wavefunctions are de¬ 
termined by the Berry connection A and Berry curvature 

B, 


3 

hy = Mis sin(rj • k^), (S.21) 

1=1 


hz = + 2ih'^n(^1+1)1 sm[rj ■ (k,+i -k,)], (S.22) 

1 = 1 

and the pseudo spin aj = (crj, aj, cr|) with the Pauli 
matrices for the jth atom defined as aj = | ) (mr^. | + 

|^r,-)(er,-|, cTj = -i\er.){mrj\ +i\mr.){er-\ and cr| = 
|er/)(erj - 




(S.23) 


B = V X A, 


(S.24) 


where V and Vx are the gradient and curl operators 
with respect to r. 

Similar to a spin-1/2 in a magnetic field h, the eigenen- 
ergies of the SL eigenstates of the effective Hamiltonian 
in equation ( S.19| ) are 


€± = 


(S.25) 





























where ± are for the upper and lower bands and h is the 
magnitude of h. The eigen wavefunction in the upper 
band can be written as 


l^+W) 


__ / h + hz \ 

^/2h(hTT^ yhx+ihy ]' 


(S.26) 


This wavefunction is well defined except for the south 
pole of the Bloch sphere where h = (0,0, with = 
—h. The eigen wavefunction near the south pole can be 
written as 


W(r))=e-.W,W) = ^=^ 

where the gauge transformation 


hx ihy \ 

h-h, 

(S.27) 


gi0(r) _ 'i'hy 

I ihy I 


(S.28) 


This way, the whole Bloch sphere is fully covered by two 
gauges [6]. 

The Berry connection in the upper band (except for 
the south pole of the Bloch sphere) is 


A+ =i(^+|V|V’+) 

1 


2h(h T 

and near the south pole we can use 


{hx^ hy hy^ hx^ , 


(S.29) 


AV=A+-V0(r). 

The Berry curvature is 

B+ =V X A+ 

= - 2zlm{dxi^^\dyi;^) 

2 ; 

— 2^3 ^^^^^^^xhbdyhc^ 

where Cabc with a^b^c G x, y, z is the Levi-Civita symbol. 
The Chern number is defined as the total Berry curvature 
in the whole first Brillouin zone 


(S.30) 


(S.31) 



dS 


1 

47r 



^3 ^abcha^xhbdyhQ dx dy^ 


(S.32) 


which counts the winding number of the effective mag¬ 
netic field h wrapping around the Bloch sphere in the 
whole Brillouin zone Hi]. 

In Fig |S.3[ we plot h in the Brillouin zone for a topo¬ 
logical nontrivial SL where the band gap is opened by the 
NNN hopping. If we go from the K point r+ = — ^^^ y 

to the K' point r_ = h moves from the north 

pole to the south pole of the Bloch sphere. The Chern 
number is one. If the band gap is opened by the on-site 






K K' 


(a) 


(b) 


Figure S.4: The eigenstates at K and K' points for topolog¬ 
ical and trivial SL’s. (a) For a topological SL, the eigenstate 
|e) at r+ (i.e., K point) is in the upper band while it is in 
the lower band at r_ (i.e., K' point), (b) For a trivial SL, |e) 
is the eigenstate in the upper band at both r+ and r_. Red 
(Blue) lines denote the bands with more components in |e) 


offset, Ac, hz is a constant in the Brillouin zone, h can 
only cover a patch on the Bloch sphere and the Chern 
number is zero. 

The two bands of iLeff have the smallest gap at the K 
and K' points r±, where 

3 

y^gir±.k, ^ 0 , (S.33) 

1=1 

and hence hx = hy = 0. At these symmetry points of the 
Brillouin zone the Hamiltonian is diagonal, 

H(r±) = ± 2V3n'na]az, (S.34) 


where a = - Yd=i sin((/>z+i - 

Let us consider the specific case that Ac = 0, = 

(/ — l)47r/3, and thus a = 3^/3/2. In this case 'H{r±) = 
:^9hQ'az has opposite signs at r±, as shown in Fig |S.3| 
(a). It is obvious that at r+, the eigenstate in the upper- 
band with eigenenergy 9hfl' is |e), while at r_, |e) is the 
eigenstate in the lower band, as shown in Fig |S.4| (a). 

Near r+, we can use equation ( S.26| ) to describe the 
wavefunction. However, equation (S.26| cannot describe 


all the wavefunctions in the upper band. It has a sin¬ 
gularity at r_ where the magnetic field h points to the 
south pole. We can remove the singularity by tuning Ac. 
When Ac > the effective magnetic field h does 

not experience the whole Bloch sphere [see Figl^(b)]. 
In this case, hz in equation (S.34) has the same sign for 
r±, |e) is in the upper band at both r±, as shown in 
Fig |S.4| (b). The wavefunction can be described by the 
same gauge. Thus the single-valued Berry connection A 
has no singularities on a closed surface and the Chern 
number is zero. Generally, the Chern number of the up¬ 
per band can be written as 

C = i sign ^Ac + — sign ^Ac — 4\/3f2W^ . 

(S.35) 

Specifically, when Ac = 0, C = sign((a). 
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DYNAMIC EVOLUTION 

We can solve the dynamics of the atomic ensemble in 
real space with the total Hamiltonian including the probe 
field, 


N 

i=i 

+ [{hx - ihy)\er.){mr-\ + h.c.] 

— \ eYj){G\ + h.c.], 

(S.36) 

with Vp = Ac /2 — Ap being the probe detuning with 
respect to the middle of the band gap and Vtp the probe 
field Rabi frequency. The wavefunction in real space 


For uniformly distributed atoms, in the limit V ^ U^/ye, 
we can assume all points in the Brillouin zone are oc¬ 
cupied by atoms, so the summation can be written as 
integration in the first Brillouin zone. 




^^■^dxdy, 


where S is the area of the first Brillouin zone. 


(S.44) 


DIFFRACTION CONTRAST AND TOPOLOGY 


In this section, we quantitatively compare the contrast 
between |ck^P and |ck_ P in topological and trivial SL’s. 
We assume z>p = 0 and 7 ^ = 0. Then Eq.(S.41) becomes 


N 

1^') = Cg\G) +y]ce(rj)|er^) +Cm{Yj)\mr^). (S.37) 

j=l 

The dynamic equations of the probability amplitudes are 


= ,2 \ (S.45) 

We assume O', Ac <C O 5 and thus Ce(r) is highly centred 
at K and K' points, where we have hx = hy = 0 and 


Ce{rj) = 


hVp) - 7e 


Ce{rj) 


(S.38) 


^{hx ^hy^Cm(x^ 


- -{hx + ihy)ce{rj). 


(S.39) 


where 'je/m is the decoherence rate of |e/m) states. 

In the limit of weak probe field Qp <C 7 e, we have 
cg ^ 1. In the steady state, Ce(rj) = Cm(rj) = 0. From 


Eq.(S.39), we obtain 


Cm(rj) = 


hx “h ihy 


hz hi^p “I” ih'yxri 


-Ce(rj). 


Substituting Eq.(|S.40|) in Eq.(|S.38|), we obtain 

CeiVj) = 




hz — hVp — ihje H- 


(S.40) 


(S.41) 


as 


hz+hUp+ihjrn 

In the SL coordinates, the wavefunction can be written 
|4^) = cg\G) + ^Cklck) + Ck-ki l^k-ki)- (S.42) 

k 

where the probability amplitude 

N 


Ck = (ekl^') = 

V ^ j = l 


(S.43) 


C ('r'l =_ ^ 

^ hz-ih^e 


(S.46) 


According to Eq.(S.44), the probability amplitude Ck can 
be approximately calculated by the integration of Ce(r) 
in small areas near r±. 


^ ^/NhQ,p ^^^(k„-k)-r+ 


{hz — ihje) ^dxdy 

(S.47) 


ghkp k)-r_ JJ 1 


where the phase factors have been taken out of the in¬ 
tegration since they do not change appreciably in those 
small areas. 

Eor a topological SL where U' 7 ^ 0 and Ac = 0, 
hz has opposite signs at r±. We denote ffj.^(hz — 
ihje)~^ dxdy = p + iq^ then ff^ {hz — ihje)~^ dxdy = 
—p + iq. Thus 

Ck {psm[{kp - k) ■ r+] + gcos[(kj, - k) ■ r+]}. 

(S.48) 

Eor kp = —ki, k+ = —ks and k_ = —k 2 , we have 
|ck+r + 7k_r 

Eor a trivial SL where U' = 0 and Ac 7 ^ 0, hz has 
the same sign at t±. We have ffj.^{hz — ihje)~^ dxdy = 

iJr {hz — 'ih^e)~^ dx dy = p -\- iq and 

Ck « 2 ^ cos[(kp - k) • r+]. (S.50) 
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Therefore, Ck+ = Ck_ and 77 = 0 . 

For a trivial SL with 7^ 0 and 7^ 0 , it is difficult to 
analytically calculate the results. Nevertheless, as shown 
in Fig .3 (b) in the main text, 77 reduces to a small value 
in the (7 = 0 phases especially for large O'. 


COUPLED- WAVE EQUATIONS 


In this section, we calculate the propagation of a 
weak probe pulse. We set = —ki, k+ = —ks and 
k_ = —k 2 . There will be directional emission in k± by 
TDS |ek±). The emitted photons in k± interact with 
the atomic ensemble the same way as the probe field k^, 
creating TDS and resulting in directional emission in the 
other two modes. For example, emission at k+ excites 
the atom ensemble and leads to emission at k_ and k^. 
Therefore, the three optical fields in k^ ± are coupled 
via the atoms. In the following, we derive the coupling 
equations. 

The expectation value of the dipole moment of the jth 
atom at position Vj is 


(^lAil^') =y]ck(G|/ij|ek)+C.C. 
k 


Vn 




+ C.C., 


(S. 51 ) 


where jlj is the dipole operator of the jth atom and 
/i = {gj\fij\ej). Since the atoms are homogeneously dis¬ 
tributed, the polarization density as a function of posi¬ 
tions is. 


P(r) 


N li 

V TV 




+ C.C., 


(S. 52 ) 



Figure S.5: The contrast of the NNN coupling strengths. 
The contrast (|fli 3 | — |fli 2 |)/(|fli 3 | + |f^i 2 |) where Q 13 (fli 2 ) 
is the NNN hopping coefficient from |ekp) to |ek^) (|ek_)). 
The dot lines are zero points. 


m i = z/p/c = kp. The coupled-wave Maxwell equations 
for the three modes are 

f ^ 1 /I \ j k 

• V + -^ j ^ iXoEp + X-E+ + X+E-), 

^ (S. 53 ) 

^ ^ I) f ■ 

(5. 54 ) 

(5. 55 ) 

where kp/^/_ are the unit vectors in the directions of 
hpij^l_ and Vg is the group velocity of the pulses. 


where V is the volume of the atomic ensemble. The 
probe field Rabi frequency is flp = jj^Ep/h with Ep be¬ 
ing the probe field strength. The polarization density 
contains the Fourier components in k±, P±^p = eox±Pp 
with x± — iic\^^/Ve{)Ep and cq is the permittivity 
in vacuum. The notation P±,p means the k± Fourier 
component of polarization density generated by optical 
field in k^. Once modes k± are excited and the corre¬ 
sponding fields E± are generated, they also polarize the 
atoms. For (ji =0, 02 = 47r/3 and 03 = 27r/3, the three 
fields have a cyclic relation between each other. We have 
= eoX-^+, P-,+ = eoX+^+, Pp- = eoX+E- and 

P+- = eoX-^-* 

We assume the quasi-static approximation in which the 
atoms are assumed in steady sates, i.e., Ck is given by 
equation (S.44). This approximation is justified when the 
pulse duration is much longer than the decoherence time 
l/Xe- We also use the slowly-varying-envelop approxi¬ 
mation. The fields in the three relevant modes kj with 
j = p, + and — are denoted as where 


LOCAL AND GLOBAL EFFECTS 

In the main text, we use the chiral topological current 
to explain the consistence between 77 and the topological 
order. The topological current is a global effect for which 
we must treat the lattice as a whole. Since the three 
superradiant TDS |ek ±) are next nearest neighbours in 
the lattice, the local effect induced by the transitions 
inside a unit cell is also important. The direct NNN 
transition strength from jck^) to |ek+) or |ek_) is Q13 
or Qi 2, respectively. If the dynamics is dominated by 
these direct transitions, we expect the contrast of the 
probability 77 should be similar to (|f^i31 — |1)/(|^131 + 
IQ12I), which is plotted in Fig ]S. 5 [ 

We plot the contrast 77 in Fig |S.6| for different 
The larger Q' is, the closer 77 is to Fig |S. 5 | and the more 
important is the local effect, as shown in Fig]^ (a) and 
(b). For O' < 7^, the effect of the direct NNN transitions 
is suppressed by the decoherence rate 7^. The global 
topological effect dominates for most areas and 77 is closer 
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Figure S.6: The global and local effects. Diffraction contrast 
T] for an SL with (a) Q' = 1, (b) Q' = 0.01, (c) Q' = 0.001 
and (d) Q' = 0.0001. = 0. The dot lines are zero points. 

Qs = 3. Ap = Ac = 0. 7e = 1 and 7m = 0.001. 


to Fig.2 (a) in the main text, as shown Figj^ (c) and 

(d). 
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